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The acoustic cloaking theory of Norris 1 permits considerable freedom in choosing the transformation 
/ from physical to virtual space. The standard process for defining cloak materials is to first define 
/ and then evaluate whether the materials are practically realizable. In this paper, this process is 
inverted by defining desirable material properties and then deriving the appropriate transformations 
which guarantee the cloaking effect. Transformations are derived which result in acoustic cloaks 
with special properties such as 1) constant density 2) constant radial stiffness 3) constant tangential 
stiffness 4) power-law density 5) power-law radial stiffness 6) power-law tangential stiffness. 7) 
minimal elastic anisotropy. 



PACS numbers: 43.20Mv, 43.40. Sk, 43.30.Wi, 43. 20. El 



I. INTRODUCTION 

Acoustic cloaking refers to making an object invisi- 
ble to sound waves. This is achieved by enclosing the 
object of interest with an acoustic cloak which guides 
waves around the object. The cloak leaves the wave- field 
outside the cloak indistinguishable from the wave-field 
without the object present. The phenomenon of cloaking 
is not restricted to acoustics but can occur for different 
types of waves such as electromagnetic waves^, elastic 
waves^, and in a more exotic example, quantum mechan- 
ical systems^. We restrict our attention here to acous- 
tic cloaking, and specifically pentamode acoustic cloak- 
ing for which the density is isotropic. The reader is re- 
ferred to the review articles by Bryan and Leise^ and 
Greenleaf et al. 6 for a comprehensive review of different 
types of cloaking, its historical development and relation 
to previous work in inverse problems^. A review dedi- 
cated to acoustic cloaking and transformation acoustics 
is provided by Chen and CharP. 

Initial work in acoustic cloakin^P^ was based on 
transformation optics as developed by Pendry et a/P^. 
Cummer et al. 9 mapped the 2D acoustic equations in 
a fluid to the s ing le polarization Maxwell's equations, 
while Chen et aZp2 mapped the 3D acoustic equation to 
the direct current conductivity equation in 3D. Cum- 
mer et al. m derived a formulation for 3D acoustic cloak- 
ing starting from scattering theory. These formulations 
achieved acoustic cloaking using anisotropic density and 
isotropic stiffness. Norris^ provided a formulation of 
acoustic cloaking which using both anisotropic inertia 
and stiffness, and as a special case, derived a formula- 
tion using isotropic density and anisotropic stiffness. 

The acoustic cloaking theory of Norris^ involves map- 
ping the physical space to the virtual space using the 
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transformation / as illustrated in Figure [T] The mate- 
rial properties of the cloak can be obtained by choosing 
the transformation / and using eqs. ^ to compute the 
material properties. In practice, however, the properties 
obtained may not be useful because they are unattain- 
able. In this paper, we derive special forms of / which 
may result in physically realizable cloaking metamate- 
rials, which are composite materials whose macroscopic 
acoustic properties are controlled by engineering their mi- 
crostructure. The design and fabrication of such acous- 
tic metamaterials is possible because of recent advances 
in material science and engineering. Cloaking metama- 
terials may have spatially varying anisotropic density 
and stiffness. We restrict our attention here to spa- 
tially varying material properties with isotropic density 
and anisotropic stiffness because Norris 1 showed that 
anisotropic density implies that the acoustic cloak has in- 
finite mass. He presented an alternative acoustic cloaking 
formulation involving pentamode materials which have 
isotropic density and a special type of anisotropic stiff- 
ness. Since cloaking is achieved with anisotropic stiffness 
as opposed to density, it is expected to have frequency in- 
dependent behavior in theory. In practice, however, the 
behavior is expected to be only extremely wideband or 
weakly frequency dependent because of frequency limita- 
tions arising from: 

• the size of the virtual cloak radius S as compared 
to the wave-length of the incident acoustic wave A. 
Cloaking is ineffective for incident acoustic waves 
whose wavelength is of the same order as the virtual 
cloak radius S. 

• the length scale of periodic structures present in 
the composite material used to fabricate the pen- 
tamode material. 

• an intrinsic frequency dependence in the properties 
of the composite material. 

The paper is organized as follows. Section [TT| pro- 
vides a short review of the pentamode acoustic cloak- 
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TABLE I. Feasibility of constant material properties in 2D 
(d = 2) and 3D (d = 3). 
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FIG. 1. The transformation / from physical space (right) 
to virtual space (left). The particular mapping shown is the 
linear mapping of eq. |2^i). 



ing theory^. Next, in ^ III , we derive transformations / 
which yield specialized spatial distributions of material 
properties, namely, 1) constant density 2) constant ra- 
dial stiffness 3) constant tangential stiffness and explain 
the wave-propagation with ray-tracing. Such distribu- 
tions may be simple to manufacture and may also help 
in evaluating the feasibility of manufacturing material 
properties on Ashby charts, as irffl. We note that in re- 
lated work, CummeP' has derived transformations for 
electromagnetic cloaking which yield constant magnetic 
permittivity e^. In j ]IV| we derive transformations which 
yield 1) power-law density 2) power-law radial stiffness 
3) power-law tangential stiffness. In we derive a dis- 
tribution of elastic properties that minimizes the elastic 
anisotropy. 



II. REVIEW OF ACOUSTIC CLOAKING USING 
PENTAMODE MATERIALS 

Acoustic cloaking relies on a transformation from an 
undeformed or original domain Q to a current (deformed) 
domain uj which is given by the point-wise deformation 
X G ft —> x G uo. Using notation from the theory of finite 
elasticity, the deformation gradient is denned F = Vxx, 
or in component form Fu = dxijdXj. The Jacobian of 
the deformation is J = deti 7 ", or in terms of volume ele- 
ments in the two configurations, J = d v/ d V. The polar 
decomposition is F — Vi2, where R is proper orthogonal 
(RR l = R f R — J, det R — 1) and the left stretch tensor 
V is the positive definite solution of V 2 = B where B is 
the left Cauchy-Green or Finger tensor B — FF l '. 

For a given transformation the cloaking material is not 
unique 1 . For instance, the inertial cloak!2(eq. (2.8)) is de- 
fined by the density tensor p = po^B -1 and bulk mod- 
ulus K = KqJ. At the other end of the spectrum of 
possible materials is the pentamode cloak with isotropic 
density, which can be chosen if the deformation satis- 
fies the property^ (Lemma 4.3) that there is a function 
h(x) for which div hV = 0. This is the case for radially 



symmetric deformations in 2D and 3D, the cylinder and 
sphere, respectively. The pentamode material is then 1 
(eq. (4.8)) p = poJ' 1 , K = K J, S = J _1 V, where the 
fourth order elasticity tensor is C = KS (8) S. 

Radially symmetric deformations in 2D and 3D are 
defined by R~ X X = r^x where R = |X|, r = |x|. If 
we let R = /(r), then the inverse mapping is defined as 
X = /(r)r _1 x. In this case we can identify a "radial 
bulk modulus" and an orthogonal bulk modulus 

K r = C rrrr = KoJ~ 1 B rr , Kj_ = K§J~ X ~Bqq, (1) 

where denotes a direction orthogonal to the radial di- 
rection. 

The requirements on the transformation generator f(r) 
admit an infinity of functions. In this paper, we will take 
advantage of this fact to design transformations which re- 
sult in desirable material properties. One familiar trans- 
formation, used frequently in acoustic cloaking work, is 
shown in Figure [l] Some examples from the infinite fam- 
ily of permissible transformations are: 
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where 6, a and S are the outer radius of the cloak, the 
inner radius of the cloak, and the virtual cloak radius, 
as shown in Figure [l] The linear transformation is 
known as the KSVW mapping for the paper 15 in which 
it was first used extensively in this form. Figure |2(a)| 
shows rays passing through a KSVW cloak. The power 
law mapping (J^Jd) yields constant p r for the inertial cloak 
and constant K r for the pentamode cloak in 2D. As we 
will see below, the transformation (J^) yields constant 
K for the inertial cloak or constant p for the pentamode 
cloak. 



III. TRANSFORMATIONS YIELDING CONSTANT 
MATERIAL PROPERTY DISTRIBUTIONS 

In this section we determine the transformations / that 
yield constant spatial distributions of material proper- 
ties. We consider a constant distribution of p in §III.A[ 
constant K r in ^ III.B and constant K± in ^ III.Q re- 
spectively. Both 2D and 3D cases are considered. The 
conditions for feasibility are summarized in Table [T[ 
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The density p (isotropic), the radial stiffness K ri and 
the tangential stiffness K± for a pentamode cloak sur- 
rounded by an anisotropic fluid with density po and bulk 
modulus Kq satisfy the following relations: 



Kr{r)=K ^(t) d -\ K ± (r) = K f^) 



d-3 
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(3) 



Our procedure to determine / consists of treating eqs. 
([3| as differential equations for / with the material prop- 
erties (p, K r ,K±) known. Having determined / we will 
prove that it satisfies the necessary conditions / > 0, /' > 
for r G [a, 6], the existence of b such that f(b) = 6, and 
the existence of 8 such that f(a) = S«a. 

We remark that eqs. (|3| are consistent with the connec- 
tion between the three parameters which is independent 
of the transformation: 



K^ d K r K d ~ x 



(p/po) d - 



(4) 



Equations <(3|) also imply that at the edge of the cloak, 
the cloak is impedance matched in the radial direction 
but not in the tangential direction: 



Z r {r) = ^K^> Z r (b) = ^/K^~o = Z , 
Z±(r) = ^KZp => Z ± (b) = Z f(b). 



(5) 



In contrast, at the edge of the cloak, the wave speeds are 
matched in the tangential direction, but not in the radial 
direction: 



c r (r) = VX/P => c r (b) = co//'(6), 
c±(r) y/K±/p => c±(b) = c . 



(6) 



A. Transformations yielding constant cloak density p 

We assume that we are given the cloak geometry b > 
a > 0, and using equation ([3| and p(r) = p c , a constant, 
we determine / and prove that / > 0, /' > for r G 
(a, b). Solving (|3| for / yields 



/= \b d + ^(r d -b d )] 
Po 



(7) 



and differentiating implies /' > for p c > 0, showing 
that / is monotonically increasing. Since f(a) = 8 > 0, 
/ > G (a, b). Enforcing 8 = f(a) < a yields p c > p Q . 
This result makes physical sense because the deformation 
/ _1 compresses the volume of fluid into a smaller volume. 
K r , K± are determined by using / in equation (|3| and 
are given by: 
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Similarly 5 and its sensitivity can be determined by 



(8) 



Finally, note that, K r ex J 2 ^" 1 ) and K± ex S~ 2 In 3D 
this implies a very strong decrease in K r with 8. 

We note that rays in the cloak are straight lines in 
deformed space 1 . This allows us to trace the rays for 
the constant density cloak by deforming the straight rays 
corresponding to a plane wave traveling through a ho- 
mogeneous medium by the inverse tra nsform ation 
Ray-tracing results are shown in Figure |2(b)| from which 
we can see that the rays curve gently in the outer region 
of the cloak and sharply close to the inner radius. In our 
experience, the smooth nature of this propagation makes 
it easy for this to be simulated with standard linear finite 
elements 16 . 



B. Transformations yielding constant radial stiffness K r 

We formulate the problem by assuming that we are 
given a > 8 > and K r constant. As before, start with 
the expression for K r from equation (|3|, and treat it as 
an ordinary differential equation for / and get 



f f, = K f dr 

J f*-i ar kJ H-i- 



(10) 



The solution is different for 2D and 3D and is determined 
separately in the next two subsections. 

1. Transformation yielding constant K r in 2D 

The transformation / in 2D can be determined by solv- 
ing equation ( 10 ) to get 



The cloak parameters then follow from eq. ^ as 
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(12) 



We prove that given a > 8 > 0, it is possible to find b 
such that b>a>8>0\/K r = constant > 0, K r ^Kn. 
Since f < 1 and a < 6, the exponent — 1 in eq. ( 12) 
must be positive to ensure b > a. i.e. < K r < Aq 
must hold for constant K r cloaks in 2D. This shows that 
the cloak outer radius is greater than the inner radius, 
making the cloak physically realistic. 

This case is of interest because the only parameter that 
varies with r is the density. At the outer radius p(b) = 
po^ 2 - > Po? and the value at the inner radius is p(a) = 
p{b){5/af. 

The rays for a constant K r cloak in 2D are shown in 
Figure 2(c)| Unlike the rays for the constant density 
cloak shown in Figure [2(b)| the rays for constant stiff- 
ness curve sharply at the outer surface of the cloak. In 
our experience, the propagation of waves near the surface 
is extremely hard to capture with standard linear, time- 
domain, finite elements, possibly due to sharp change in 
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Ray tracing for KSVW Mapping 
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(a) KSVW cloak 



Ray tracing for Constant Rho Mapping 




(a) Initial propagation. 



(b) Constant density cloak 



(b) Complete propagation. 



Ray tracing for Constant Kr Mapping 
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(c) Constant stiffness cloak 



Ray tracing for Optimal mapping 




(d) Cloak with optimal anisotropy 

FIG. 2. Rays for KSVW, constant density, constant stiffness, 
and optimal anisotropy cloaks in 2D, a = 1, b = 3. 



FIG. 3. Rays for power law K r , with parameters a = 1, b = 3, 
a = 3. 



the direction of propagation. The required element den- 
sity is in the order of hundreds of elements per shortest 
wavelength. This is contrast to the rule of thumb in tran- 
sient, explicit finite element analysis in which typically 
sixteen elements per shortest wavelength are used. 



2. Transformation yielding constant K r in 3D 



Treating equation (10) with d = 3 as a differential 



equation for / with K r known and constant, gives 
-1 Kn/1 In-i- 1 



= f ({) 2 >°- ( 13 ) 



J l b K r \r b) 
Using the above definition of / in eq. ([3| we get 
K ( bK r 



P = Po- 



K r \bK + r(K r -K ) 



K± = K r 



bK 



5 = /(a) = 



bK + r(K r - K ) 
abK r 



bK + a(K r -K )' 



(14a) 

(14b) 
(14c) 



We now prove that under certain conditions it is possible 
to find b such that b > a > S > 0, meaning that the cloak 
outer radius is greater than the cloak inner radius. As a 
consequence we prove that / > 0. We start by rewriting 
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1 2 3 




(a) Initial propagation. 



(a) Initial propagation I. 





(b) Complete propagation. 



(b) Initial propagation II. 



FIG. 4. Rays for power law p, with parameters a = 1, b = 3, 
a = 3. 



equation (14c) as 



(15) 



To ensure | > 1, we require if r < if of • Since, /' > G 
(a, 6) and /(a) = (5 > 0, it follows that / is monotonic 
in the interval (a, 6) and therefore / > G (a, 6). Note 
that, if 5 « a ==> if r << if • Since i^i = K$/K r , 
this means >> if r . Thus, to ensure a small scatter- 
ing cross section we need a material which is very stiff in 
one direction (_L) as compared to the other (r). We ex- 
pect that such a material will require careful engineering, 
and such practical constraints may restrict the amount 
of reduction in the scattering cross section that can be 
achieved. 



C. Transformations with constant tangential stiffness K± 

1. Transformations yielding constant K± in 2D 

In two dimensions, constant K± is the same as con- 
stant K r (which was previously considered), because 
K r K ± = Kq in 2D. 




(c) Complete propagation. 

FIG. 5. Rays for power law />, with parameters a = 1, 
a — —3. 



2. Transformations yielding constant tangential stiffness in 3D 

We formulate the problem as follows. We consider that 
we are given a > S > 0. We need to find a / > 0, /' > 
satisfying f(a) = 5, f(b) — 6, b > a. We treat eq. Q as 
a differential equation for / and obtain after using the 
boundary condition f(a) = J, 



f(r)=5+^(r-a)>0 and/ 



K± 



> 0, (16) 



for r > a. Using b = f(b) > a yields the condition 
K± > Kq. The constant K± cloak in three dimensions 



5 



ations of density, radial stiffness tangential stiffness and 
a case in which density is proportional to stiffness. This 
treatment is restricted to two dimensions. 
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(a) Initial propagation. 




(b) Rays converging. 
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(c) Complete propagation. 
FIG. 6. Rays for an acoustic concentrator f(r) — b(r /b) 1 ^ 10 . 



is characterized by / of eq. (16) and 



K ± ,5 K ± y 



(17) 



Note that this is the same transformation as the KSVW 
transformation in eq. . 



IV. TRANSFORMATIONS YIELDING POWER LAW 
PROPERTY VARIATION IN 2D 

We now consider more complicated spatial distribu- 
tions for the material properties, namely power law vari- 



A. Transformations yielding power law density 

Consider the following strategy: Given b > a > 5 > 0, 
a ^ 0, we determine p(a) = p a that ensures f(b) = b 
and prove that p a > 0. The power law for density we 
consider is as follows: 



p(r) 



f'dr 
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-dr, (18) 



where the latter is a consequence of eq. (J3|. The cases 
a ^ — 2 and a = — 2 need to be considered separately. 
To summarize, the power law density cloak in two dimen- 
sions is characterized by (18)i, with 
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(21) 



Note that the a = corresponds to the special case of 
constant density in 2D, considered in §111. A| To prove 
p a > for a ^ — 2, consider the two cases a > — 2 and 
a < —2. In the first case we have, b a+2 > a a+2 , and 
therefore p a > 0. In the second case, we have, b a + 2 < 
a a+2 . In addition, a+2 < and therefore (a+2)/(6 a+2 - 
a a+2 ) > 0. Hence p a > in this case as well. The 
same argument, substituting r for b can be used to prove 
f 2 > 0, and choosing the positive root, we get / > 0. 
The positivity of /' follows. A similar analysis can be 
performed for a = —2. 



B. Transformations yielding power law radial stiffness K r 

Here, we consider the following power law for K r 

K r (r) = K a {£j° , a/0. (22) 



The case a = (constant K r ) was considered in j |III.B| 
Given 6>a>(5,a^O, we find what value of K a ensures 
that f(b) = b. We then prove this K a > 0. Using eqs. (|3|, 
(22) plus the boundary conditions f(a) = S and f(b) = 6, 
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we get 

-(;)")}• 

if o" 

if r a+l' 

-(?)*)• 

/ 2 ifo a a 
each of which is clearly positive. 

C. An acoustic concentrator 



(23a) 
(23b) 
(23c) 

(23d) 



Motivated by the form of equation (23a), we con- 
sider the transformation r = b(R/b) 10 or equivalently, 
R — f( r ) = ^(r/6) 1 / 10 , yielding the rays shown in Figure 
[6] The focus of the rays can be made arbitrarily tight. 
The process can also be reversed: one can place a source 
at the focus and convert a cylindrical wavefront gener- 
ated by a point source into a plane wavefront. Similar 
work on designing acoustic concentrators using transfor- 
mation acoustics has been recently reported in Wang et 
a/P^l Previously, Rahm et al 18 reported the design of 
electromagnetic concentrators. 



D. Transformations yielding power law tangential stiffness 



Requiring that < S < a implies the a range of possible 
values for j3. The special case of A = 0, i.e. a = (d— 2)/3, 
needs to be distinguished. In 2D, this means a = and 
therefore K r = constant, and is therefore not interesting. 
In 3D, this implies, a = ~, and eq. p5| can be integrated 
to give a power law solution for /(r)T 



t-H 



[I 



(d = 3, <*=-), (29) 



In this case, S = b (§) M , which is clearly positive, satisfies 
the constraint S < a only if \i > 1, or equivalently, j3 < 

1/3 

Kq/pq , setting an upper limit on f3. 



In summary, eq. ( 24 ) has cloak-like solutions for a > 
1 in 3D, with associated limits on the 



in 2D, and a > 3 
possible range in value of the parameter j3. 



V. TRANSFORMATIONS YIELDING MINIMAL ELASTIC 
ANISOTROPY 

Here, we consider acoustic cloaks which have minimal 
elastic anisotropy in a certain sense. We are motivated 
by the fact that extremely anisotropic materials are hard 
to design and manufacture. Minimizing anisotropy there- 
fore may lead to a practical cloak. We begin by defining 
two measures of anisotropy in equation (31) and prove 
that only one of them yields physically meaningful trans- 
formations. 



In two dimensions, this is equivalent to power law K r , 
because K± = Kq/K v . 



E. Transformations yielding proportional density and radial 
stiffness 

Since density is usually associated with stiffness, we 
consider a power law linking density with stiffness in the 
radial direction. This power law is defined as 



(24) 



for a, /?(> 0) constant. Using the pentamode relations 
([3]) for K r and p in the above equation, we get a differ- 
ential equation for / 



f^f =u x r x ~ 1 



where 
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When A ^ eq. (25) yields 



f(r) = (V+u; A (r A -a A )) 



l/A 



(25) 
(26) 

(27) 



This equation is not consistent with the constraint that 
f(b) = b > a if A < 0, in general. For A > 0, setting 
f(b) = b yields 



(b x -L0 X (b x -a x )y 



(28) 



A. Optimal transformations in cylindrical cloaks 

We define the following parameter to be a local mea- 
sure of the anisotropy in the cloak, 



7 : 



a 



1 
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where a 



(30) 



This is the same anisotropy parameter introduced irP. 
It can be shown that the minimum value of eq. ( 30 ) is 2 



and it occurs for a = 1, i.e. when there is no anisotropy. 
Based on eq. (30) we introduce two global measures of 
cloak anisotropy, 



7i 



= VJ 1 [ 7 dF, 



72 



7 dv, 



(31) 



where V u , Vq, are the volumes (areas) in the physical 
and virtual domains, respectively. It follows from eq. ([I]) 
and the identity Q for d=2 (K$ = K r K±) that 



J _1 trB = 
Similarly, 

JtrB" 1 = 



K r 
K~ 



K± 
K 



K r 



K± 

K r 



7- 



(32) 




7- (33) 



Based on the identities (32) and ([33]), it follows that 

71 = KT 1 L trB" 1 dv = Vj 1 LfJfJ dv, (34a) 

72 =V^ / n trB dV = V^ J n §^§^ dV. (34b) 
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The parameter 71 is therefore the average in the current 
configuration of the sum of the principal stretches of the 
mapping from the original (virtual) domain. Conversely, 
72 is the average in the original configuration of the sum 
of the principal stretches of the inverse mapping from the 
current (spatial) domain. 

The global anisotropy measures 71 and 72 are min- 
imized by the Euler-Lagrange equations. Consider 71, 
then assuming lj is fixed, we have 



<5 7 i = 2V- 1 j 



dv 



2 



duo 



dXi 5( dXj 

dxj ^ dxj 

— — oXi rijds 
dxj 



d 2 X z 
dxjdxj 



^jdv. 



The surface integral vanishes because, by assumption, 
the value of X on the boundary of is constant (in fact 
X = x is required on dui), and therefore we deduce 



mhi7i <^> V 2 X = in lj, 
min72 <^> V^x = in Q. 



(35a) 
(35b) 



It is interesting to note that these equations are satisfied 
by conformal transformations, a large class of potential 
transformations. Here, however, we restrict attention to 
purely radial transformations. 

Consider ( |35a| ) first. Assuming the inverse mapping 
X = f(r)r~ x then it is straightfor ward to show that 



V 2 X = [r(rf')' - fjifr 2 )- 1 ^ and ([35a]) is satisfied if 
/ = At + Br~ x , for constants A and B. As before, 
we assume the cloak occupies R G [5,6], r G [a, b] with 
< S < a < b. The constants are then found from the 
conditions f(a) = 5 and f(b) = b, yielding 



/ (r) = (b 2 - a 2 )" 1 



(b 2 - aS)r - (a - S)b 2 



(36) 



The same result can be found by noting that the 



anisotropy parameter of (34a) reduces for radially sym- 
metric transformations to 



d 



7i 



(f) 2 + (d-l) 



^dr, (37) 



for d = 2 or d = 3. The minimizer satisfies the Euler- 
Lagrange equation r 3 - rf ( r rf - 1 y / ) / _ (^ _ \)f — o, which 
for d = 2 gives ( [36] ) . In the same way, we find that ( 35b ) 
is satisfied if r = AR + BR -1 , for constants A and B. 
The end conditions r(S) = a and r(b) = b imply that the 
transformation which minimizes 72 is 



r = (b 2 - s 2 y 



(b 2 - aS)R +(a- 5)b 2 



R 



(38) 



However, this transformation function is generally not 
one-to-one. The problem is illustrated in Figure [7] and 
comes from the fact that drj dR = at some value of 
R G (5,6). This cannot occur for the mapping function 
( [36] ) . Consequently equation (36) is a valid transforma- 
tion for acoustic cloaking while equation (38) is not. 



a: 0.5 




FIG. 7. The solid curve shows the transformation defined by 
(36} for {6, a, b} = {0.1,0.8,1.0}. The dashed curve is the 
mapping ([38} for the same cloak parameters. 



B. Optimal transformation for spherical cloaks 



We now take (34) as the definition of the global 
anisotropy measures. Using again the inverse map- 
ping X = /(r)r- 1 x it follows that V 2 X = [(r 2 f) f - 
2/](/r 2 ) _1 X. The transformation which minimizes 71 is 
therefore 



f(r) = (6 3 - a 3 )" 1 



(b 3 - a 2 8)r - (a - S)b 3 



The transformation which minimizes 72 is 



r = (b 6 - S 6 ) 



3\-r 



(6 3 - aS 2 )R + (a - S)b 3 



R 2 



■ (39) 



(40) 



but this again has the unphysical nature found for the 
2D case. We conclude that minimization of 72 using a 
single valued function does not appear to have a single 
or unique solution. 



C. Numerical examples 

The minimizing value of 71 may be found by integrat- 
ing (34a) by parts, and using (|35a|), 



7ln 



^ra) (2*r//')|*, 2D, 
1 



..) {^ :2 i r)': r 3D. (41) 



Thus, 



7i 




-a 2 )- 2 [(b 2 -a5) 2 + (a-5) 2 b 2 ], 2D, 
-a 3 )- 2 [(b 3 -a 2 5) 2 + (a-8) 2 b% 3D. 

(42) 

The relative value of the anisotropy parameter 71 is 
shown in Figure [8] for the three mappings of equation 
In all cases, the value of 71 exceeds the minimum 
7imin for the optimal transformations in equations (36) 



and (|39[). The KSVW mapping in 2D has anisotropy only 
slightly more than the minimum, but for 3D the KSVW 
displays much larger global anisotropy. 
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5 

(a) Two dimensions 




5 

(b) Three dimensions 



FIG. 8. The relative value of the global anisotropy param- 
eter for the three transformations denned by equation ^ 
for {a, b} = {0.8,1.0}. The curves show log(Ti — 1) where 
Ti = 7i/7imin, with 7imin given by J42| and 71 calculated 
based on the mappings in equations |2|. 



VI. CONCLUSION AND DISCUSSION 

Transformation acoustics, like its close analog trans- 
formation optics, possesses a huge freedom in the way 
that the transformation can be chosen. This paper sheds 
some light on potential choices. We have shown that is 
possible to always fix at least one of the three material 
parameters relevant to radially symmetric deformations. 
Starting from the theory of Norris^, we have derived sev- 
eral forms of the transformation / which yield specialized 
distributions of material properties such as constant and 
power law density, radial stiffness, and tangential stiff- 
ness. This was achieved by reinterpreting the governing 
equations for the material properties as differential equa- 
tions for the transformations. We derived a functional 
form of / that minimizes elastic anisotropy in a certain 
sense. 
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